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Abstract 

We study a simplified system of the original Ericksen-Leslie equations for the 
flow of nematic liquid crystals. This is a coupled non-parabolic dissipative dynamic 
system. We show the convergence of global classical solutions to single steady states 
as time goes to infinity (uniqueness of asymptotic limit) by using the Lojasiewicz- 
Simon approach. Moreover, we provide an estimate on the convergence rate. Finally, 
we discuss some possible extensions of the results to certain generalized problems 
with changing density or free-slip boundary condition. 
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1 Introduction 

We consider the following hydrodynamical model for the flow of nematic liquid crystals 
(cf. [13,15]) 

vt + v - Vv - uAv + VP = -AV ■ (Vd© Vc?), (1.1) 
V-v = 0, (1.2) 
dt + vVd = 7(Aci-/(ci)), (1.3) 

in $7 X M+, where C (n = 2, 3) is a bounded domain with smooth boundary F. Here, 
V is the velocity field of the flow and d represents the averaged macroscopic/continuum 
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molecular orientations in M" (n = 2,3). P{x,t) is a scalar function representing the 
pressure (including both the hydrostatic and the induced elastic part from the orientation 
field) . The positive constants z/, A and 7 stand for viscosity, the competition between 
kinetic energy and potential energy, and macroscopic elastic relaxation time (Debroah 
number) for the molecular orientation field. We assume that f{d) = VF{d) for some 
smooth bounded function F : R" ^ M. VdQVd denotes the nx n matrix whose (i, j)-th 
entry is given by Vjrf ■ VjC?, for I < i, j < n. 

In this paper we deal with the system (ll.ip - (11.3p subject to the initial conditions 

v\t=o = Vq{x) with V ■ fo = 0, d\t=o = do{x), for x G fi, (1.4) 

and the Dirichlet boundary conditions: 

w(x,t) = 0, d{x,t) = do{x), for (x,t) G r X M+. (1.5) 

In [13], the author proposed equations (ll.ll) - (ll.3l) as a simplified system of the original 
Ericksen-Leslie system (cf. [2,12]). By Ericksen-Leslie's hydrodynamical theory of the 
liquid crystal, the (simplified) system describing the orientation as well as the macroscopic 



motion reads as follows (here we assume the density to be constant) 

Vt + v-Vv-uAv + \/P = -AV- (Vc/0 Vrf), (1.6) 

V-v = 0, (1.7) 

d^ + vVd = -f{Ad+\\/d\'^d), \d\ = 1. (1.8) 



In order to avoid the gradient nonlinearly in (11.81) . usually one uses the Ginzburg-Landau 
approximation to relax the constraint \d\ = 1. The corresponding approximate energy is 

lll'^dl' + ^M-^fdx, 

where 77 is a positive constant. Then we arrive at the approximation system (ll.ip -l fOll . 
where 

f{d) = \{\d\'-l)d (1.9) 

with its antiderivative 

F(d) = ^(|rfp-1)^ (1.10) 

The Ericksen-Leslie system is well suited for describing many special fiows for the 
materials, especially for those with small molecules, and is wildly accepted in the engi- 
neering and mathematical communities studying liquid crystals. System (ll.ll) - (ll.3l) can 
be possibly viewed as the simplest mathematical model, which keeps the most impor- 
tant mathematical structure as well as most of the essential difficulties of the original 
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Ericksen-Leslie system (cf. [15]). System (ll.ip - (ll.3p with Dirichlet boundary conditions 
has been studied in a series of work not only theoretically (cf. [15,16]) but also numeri- 
cally (cf. [17,19]). In particular, in [15], the authors proved the existence theorem for the 
weak solutions of system (ll.ip -f |T75l ) by a modified Galerkin scheme. After that they also 
obtained the global existence and uniqueness of classical solutions to the same system for 
n = 2 or n = 3 with large viscosity assumption. Moreover, a preliminary analysis of the 
asymptotic behavior of global classical solution was also given in [15]. In the final remark 
of [15], a natural question on the uniqueness of the asymptotic limit was raised. This is 
just the main goal of the present paper. For the sake of simplicity, in the following text, 
we always treat the nonlinearity / of form (11.91) . However, it is not difficult to verify that 
our results holds true for more general nonlinearities which is analytic and with proper 
growth and dissipation assumptions. 

In this paper, we show the convergence to equilibrium of global classical solutions to 
system (ll.ll) - (|1.5p . Namely, we obtain the following results: 

Theorem 1.1. When n = 2, for any vq G Hq{^1) with V ■ f o = and do G H'^{^1), the 
unique classical solution to problem (ll.ip - (|1.5p has the following property 

lim iHt)\\Hi + \\dit)-d^\\H2)=0, (1.11) 

t— >+oo 

where d^o is a solution to the following nonlinear elliptic boundary value problem: 

-Arfoo + /(rfoo) = 0, xen, ^^^^^ 
doo = do{x), xeT. 

Moreover, there exists a positive constant C depending on Vq, d^, fi, d^o, such that 

\\v{t)\\Hi + \\d{t)~d^\\H2<C{l + t)-T^, Vt>0, (1.13) 

with 9 G (0,1/2) being the same constant as in the Lojasiewicz-Simon inequality (see 
Lemma \2l\ below). 

When the spacial dimension is three, we deal with two cases. The first result is 
concerning the large viscosity case, we have 

Theorem 1.2. When n = 3, for any vq G -f^o(f^) with V ■ vq = 0, do e H'^{Q,) under 
the large viscosity assumption v > z/o(A, 7, fo, rfo), the unique global classical solution of 
problem fll.ip - (ll.5p enjoys the same properties as in Theorem \l.l\ 

The second one is a "stability" result for the near equilibrium initial data in the three 
dimensional case. 
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Theorem 1.3. When n = 3, let d* G H'^{Vt) he an absolute minimizer of the functional 



^ Jn 

in the sense that E{d*) < E{d) whenever d = d* = do{x) on T. There is a constant a 
which may depend on A,7, and VQ,do, such that if HfoHi^i + H^o — d*\\H2 < a, then the 
problem (ll.ip - (|1.5p admits a unique global classical solution enjoying the same properties 
as in Theorem \l.l[ 

Remark 1.1. Theorem I J .31 implies that if the initial data is sufficiently close to an ab- 
solute minimizer of functional E, then there exists a global solution and the solution will 
converge to an equilibrium which may not necessarily he the original minimizer. This is 
because the set of equilibria might be a continuum. Theorem \1.3\ gives the "uniqueness" 
of asymptotic limit of the global solution to problem p.ip -( fT75D . This improves the result 
stated in [15, Theorem C], in which only sequence convergence for director field d was 
obtained. 

The problem about uniqueness of asymptotic limit for nonlinear evolution equations, 
namely whether the global solution will converge to an equilibrium as time tends to 
infinity, has attracted a lot of interests of mathematicians. If the space dimension n>2, 
it is known that the structure of the set of equilibria can be nontrivial and may form 
a continuum for certain physically reasonable nonlinearities. The reader is referred, for 
instance, to [7, Rem. 2.3.13], where the following two-dimensional equation —Am + — 
Xu = 0, \ > 0, endowed with a standard Dirichlet homogeneous boundary condition, 
is considered. And we note that, for the vector functions, the situations may be even 
more complicated. If this is the case, it is highly nontrivial to decide whether or not 
a given bounded trajectory converges to a single steady state. In 1983, L. Simon [23] 
made a breakthrough that for a semilinear parabolic equation with a nonlinearity f{x,u) 
being analytic in the unknown function u, its bounded global solution would converge 
to an equilibrium as t — oo. Simon's idea relies on a generalization of the Lojasiewicz 
inequality (see [21,22]) for analytic functions defined in finite dimensional space M™. Since 
then, his original approach has been simplified and applied to prove convergence results for 
many evolution equations (see e.g., [4-6,8-11,14,27-29] and the references cited therein). 
For our problem p.ip - p.Sp . in order to apply the Lojasiewicz-Simon approach to prove 
the convergence result, we need to introduce a suitable Lojasiewicz-Simon type inequality 
for vector functions with nonhomogeneous Dirichlet boundary condition (cf. Lemma [2TT]) . 

As far as the convergence rate is concerned, it is known that an estimate in certain 
(lower order) norm can usually be obtained directly from the Lojasiewicz-Simon approach 
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(see, e.g., [9,30]). Then, one straightforward way to get estimates in higher order norms 
is using interpolation inequalities (cf. [9]) and, consequently, the decay exponent deteri- 
orates. We shall show that by using suitable energy estimates and constructing proper 
differential inequalities, it is possible to obtain the same estimates on convergence rate in 
both higher and lower order norms. Our approach in some sense improves the previous 
results in the literature (see, for instance, [9,30]) and it can apply to many other problems 
(cf. [5,6,28,29]). 

The remaining part of this paper is organized as follows. In Section 2, we introduce 
the functional setting, some preliminary results as well as some technical lemmas. Section 
3 is devoted to the two dimensional case. We prove the convergence of global solutions 
to single steady states as time goes to infinity and obtain an estimate on convergence 
rate. In Section 4, we consider the three dimensional case. The same convergence result 
was proved for two subcases, in which the global existence of classical solutions can be 
obtained. In the final Section 5, we discuss some possible extensions of our results to 
certain generalized problems with changing density or free-slip boundary conditions. 

2 Preliminaries 

First, we introduce the function spaces we shall work on (cf. [15,26]): 

H^{n) = the closure of Co°°(n,M'') in the norm 

H-\n) = the dual of H^{n), 
H\n) = {ve L2(^],M") I v,^,v,^,^ e L2(fi,M") 

V = Co~(f^,M") n{v:V -v = 0}, 
H = the closure of V in L'^{n,W), 
V= the closure of V in H^{n), 

V = the dual of V. 

Global existence and uniqueness of classical solution to system (|l.ip - (ll.5p has been 
proven in [15, Theorem B]. More precisely, we have 

Proposition 2.1. Problem f ll.ll) -( fT75l ) admits a unique global classical solution {v,d) 
provided that Vq G i7o(f2),(io £ H'^[Q) either n = 2 or n = 3 with the large viscosity 
assumption v > z/o(A, 7, fo, (io)- 

For any classical solution {v, d) E Q x [0, T] = Qt (0 < T < +00) of problem (11. ip - 



iVfl dx 



1 < i,j < n}, 
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(II .Sp . we consider the functional 

m = limr + ^iivd^r +^ [ Pid^x. (2.1) 

It has been shown in [15] that our system p.ip - p.Sp has the following basic energy law, 
which can be viewed as a direct consequence of the balance laws of the linear momentum 
(ll.ip and angular momentum (ll.3p : 

pit) + r^\\Vv{t)r + Xj\\Ad{t) - fidm' = 0, < t < T. (2.2) 

(12 ■2p reflects the energy dissipation property of the flow of liquid crystals. Moreover, one 
can verify that £(t) serves as a Lyapunov functional for problem (ll.ip - (ll.5p . 
Next, we look at the following elliptic boundary value problem 

-Ad + f{d) = o, xen, ^ ^ 

j \ J ^ (2.3) 
d = do{x), X G r. 

Denote 

E{d) = l\\Vdf + I F{d)dx. (2.4) 

It is not difficult to see that the solution to ( 12.31) is a critical point E{d), and conversely, 
the critical point of E{d) is a solution to (12.31) (cf. [27,29] and references cited therein). 
Besides, regularity of the solution to (12. 3p has been shown in [15] such that d is smooth 
on VL provided d^ is smooth on V. 

As mentioned in Introduction, in order to apply the Lojasiewicz-Simon approach to 
prove the convergence to equilibrium, we have to introduce a suitable Lojasiewicz-Simon 
type inequality related to our present problem. In particular, we have 

Lemma 2.1. [Lojasiewicz-Simon Type Inequality] Let ip be a critical point of E{d). 
There exist constants 9 G (0, \) and /5 > depending on ip such that for any d G H^{Q) 
satisfying d\r = d^ix) and \\d — ipUn^ < l3, there holds 

\\-M + f{d)\\H-^ > \E{d)~Em'-'. (2.5) 

Remark 2.1. The above lemma can be viewed as an extended version of Simon's result [23] 
for scalar function under the use of L"^ -norm. We can refer to [10, Chapter 2, Theorem 
5.2], in which the case for vectors subject to homogeneous Dirichlet boundary condition 
was considered. Here we observe that, our present (nontrivial) boundary data for director 
field d does not depend on time. As a result, every solution to the corresponding stationary 
problem ( 12. 3p . which is a critical point of E{d) satisfies the same boundary condition as the 
solution to the evolution problem. Therefore, we only have to derive a Lojasiewicz-Simon 
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type inequality for functions d, which satisfy d\r = dQ{x) and fall into a properly small 
neighborhood of certain but arbitrary critical point of E{d). In this case, it is always true 
that the difference d = d — ip E Hq{Q). Keeping this fact in mind, we are able to prove 
the present lemma following the steps in [10, Chapter 2, Theorem 5.2] or [8]. Hence, 
the details are omitted here. We could also refer to a related case treated in [5], that a 
nonhomogeneous (time- dependent Dirichlet) boundary condition was removed by a proper 
variable transformation (cf. also [27] where the boundary condition contains a nonzero 
constant). 

In the following text, we will use the regularity result for Stokes problem (cf. [25]) 

Lemma 2.2. Denote the Stokes operator by S, which is a unbounded operator in H of 
domain H'^{Vl) n V: 

Su = -Au + VTxe H, Vm G H'^in) n V. 

Then there exists a constant C such that for any u G H'^{Q) fl V , 

\\u\\h-^ + ||7r||j/i\i; < CIISmII- 

Before ending this section, we introduce the following lemma which is useful in the 
study of large time behavior of solutions to evolution problems. We will apply it to obtain 
uniform (higher order) estimates of the solution and decay of the energy dissipations of 
system (frTD- (fL5|l . 

Lemma 2.3. [30, Lemma 6.2.1] Let T be given with < T < +oo. Suppose that y(t) 
and h(t) are nonnegative continuous functions defined on [0,T] and satisfy the following 
conditions: 

-r < Ciy^ + C2 + h{t), with / y(t)dt < C3, / h{t)dt < C4, 
di Jo Jo 

where Ci{i = 1,2,3,4) are given nonnegative constants. Then for any r G (0,T), the 
following estimates holds: 

y{t + r) < (^^ + C2r + C4) e^'^'', Vt G [0,T-r]. 

Furthermore, ifT = +00, then 

lim y{t) = 0. 
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3 Convergence to Equilibrium for Two Dimensional Case 



In this section, we prove the convergence of global solutions to single steady states as time 
tends to infinity for 2-D case. Since parameters A, 7, v do not play crucial role in the 2-D 
case, we set A = 7 = 1/ = 1 in this section for the sake of simplicity. 

When the space dimension equals to two, an important property for the global solution 
to problem (ll.ip - fll.5l) is the following high order energy law, which played a crucial role 
in the proof of global existence result in [15]. Denote 

A{t) = \^v{t)f^\\m)-im)f- (3.1) 

Then we have 

Lemma 3.1. (cj. [15, (4-9)]) In 2-D case, the following inequality holds for the classical 
solution {v,d) to problem (11.11) - fl 1.51) 

j^A{t) + iWAvf + II V(Ad - f{d))f) < C{A\t) + A{t)), V t > 0, (3.2) 

where C is a constant depending on f,^, W'^oW, \\do\\H'>-{Q)- 



3.1 Convergence to Equilibrium 

Based on the high order energy law (13.21) . we are able to show the convergence of the 
velocity field v first. 

Lemma 3.2. For any t > 0, the following uniform estimate holds 

\\vit)\\Hi + \\d{t)\\H2<C, (3.3) 
where C is a constant depending on /, ||fo||Hi) IMo||/i'2(n)- Furthermore, 

lim {\\v{t)\\H^ + II - Ad{t) + f{dm) = 0. (3.4) 

t— >+oo 

Proof. It follows from the basic energy law (12.10 that 

S{t)+ [ A{T)dT = S{0) < 00, Vt>0. (3.5) 
Jo 

By the Young inequality < \a'^ + |, we can see that £{t) is bounded from below by a 
constant which is only dependent of |n|. As a result, 
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A{t)dt < ^(0) < +00, (3.6) 
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and 

8{t)<8{Q), Vt>0. (3.7) 
fl3.7p implies the uniform estimate 



V\\d{t)\\H^<C, Vt>0. (3.8) 
Furthermore, (13. 6p together with Lemma [3TT] and Lemma [2731 yields that 

hm (||Vt;(t)|| + ||Arf(t)-/(rf(t))||)=0. (3.9) 

r— >+oo 

By the Poincare inequality, we prove the conclusion (13.41) . Concerning the uniform bound 
(13. 3p . we take r = 1 in Lemma 12731 to get 

\\Vv{t) II + II - Arf(t) + f{d{t))\\ <C, V t > 1, (3.10) 

where C does not depend on t. On the other hand, for any t G [0, 1], it follows from (13.21) 
and the fact A{t)dt < C that 



sup A{t) < e^o ^(*)'^M(0) + C <C. (3.11) 

0<t<l 

Besides, from the continuous embedding ^ Lp(1 < p < oo) and (13.81) we have 

||Ad|| < II - Arf + f{d)\\ + \\f{d)\\ <\\-Ad + f{d)\\ +C{1 + IMIlie) < C. (3.12) 

Now we can conclude (13.31 ) from (13.10I )- (I3.12I ). The proof is complete. □ 

Let S be the set 

S = {{0,u) \ -Au + f{u) = 0, in Q, u\r = do{x)}. 

The w-limit set of {vo,do) eV x H'^{Q) C L'^{Q) x if^(fi) is defined as follows: 

uj{{vo,do)) = {{voo{x),doo{x)) I there exists {t„} /" oo such that 

{v{x,tn),d{x,tn)) {voo{x) , doo{x)) iu X , as tn +oo}. 

We infer from Lemma 13.21 that 

Proposition 3.1. u!{{vo,do)) is a nonempty bounded subset in H^{Q) x H'^iVt). Besides, 
all asymptotic limiting points (foo, c^oo) of problem p.ll )-( fT75l ) belong to S. In other words, 
Lj{{vo,do)) C S. 
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In what follows, we prove the convergence for director field d. For any initial datum 
{vo,do) G V X H'^(Q)^ it follows from Lemma 1X2] that \\d\\jj2 is uniformly bounded. Since 
the embedding ^ is compact, there is an increasing unbounded sequence {t„}„gN 
and a function doo such that 

lim \\d{tn)-d^\\Hi=0. (3.13) 
In particular. Proposition 13.11 implies that cioo satisfies the equation 

- Ad^ + f{d^) =0, xen, d^\r = do. (3.14) 

We prove the convergence result following a simple argument introduced in [11], in 
which the key observation is that after a certain time to, d{t) will fall into a certain small 
neighborhood of d^o and stay there forever. 

From the basic energy law (12. ip . we can see that S(t) is decreasing on [0, oo), and it 
has a finite limit as time goes to infinity because it is bounded from below. Therefore, it 
follows from (13.131) that 

lim S{tn) = E{d^). (3.15) 

On the other hand, we can infer from (12. ip that S{t) > E{doo), for all t > 0, and the 
equal sign holds if and only if, for alH > 0, f = and d solves problem (I3.14p . 

We now consider all possibilities. 
Case 1. If there is a to > such that at this time S{to) = E{doo), then for all t > to, we 
deduce from (12. ip that 

||Vt;||=0, \\- Ad + f{d)\\=0. (3.16) 
It follows from ( II. 3p . ( I3.16p and the Sobolev embedding Theorem that for t > to 

< \\dt\\ < \\v ■ Vd\\ + \\-Ad + f{d)\\ < \\v\\L4Vd\\Li < C\\Vv\\ = 0. (3.17) 

Namely, d is independent of time for all t > to- Due to ( I3.13p . we have d{t) = doo for 
t > to. 

Case 2. For all t > 0, £{t) > E{doo)- First we assume that the following claim holds 
true. 

Proposition 3.2. There is a to > that for all t > to, \\d{t) — doo\\m < P- Namely, for 
all t > to, d(t) satisfies the condition in Lemma \2. 1[ 

In this case, it follows from Lemma ETT] that 

\E{d) - E{doo)\'-' <\\-Ad + fid)\\H-i <\\-Ad + f{d)\\, \ft>to. (3.18) 
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The fact 6 e (0, \) implies < 1 - ^ < 1, 2(1 - ^) > 1. As a consequence, 

||^;f(l-e) = ||y||2(l-'')-l||^;||<C||^;||. 
Then we infer from the basic inequality 

{a + hf-' <a^-' + h'-\ Va,6>0 

that 

{£{t) - E{d^)Y-' < (h\vf + \E{d) - E{dS 



< Q)' \\vr^'-'^ + \\-M+fm 

< C\\v\\ + \\- Ad + fid)\\. (3.19) 



Therefore, a direct calculation yields 



j^isit) - Eid^)r = -e{s{t)-E{d^)r-'j^s{t) 

^ ce{\\Vv\\ + \\-Ad + f{d)\\f 



C\\v\\ + \\-Ad + fid)\\ 
> Ci(||Vt;|| + II - Arf + /(rf)||), Vt>to, (3.20) 

where Ci is a constant depending on fo, cio, ^• 
Integrating from to to t, we get 

(Sit) - Eid^)f + C, !\\\Vv{r)\\ + II - Ad{T) + f{d{r))\\)dT 

J to 

< {S (to) - E{d^)f < oo, \/t>to. (3.21) 
Since £{t) — E{doo) > 0, we conclude that 

(||Vt;(r)|| + II - Arf(r) + /(rf(r))||)rfr < oo. (3.22) 

On the other hand, it follows from equation 01. 3p that 

\\dt\\ < ||t;-Vrf|| + ||-Arf + /(rf)<||t;|U4||Vrf|U4 + ||-Arf + /(rf)|| 

< C||Vt;|| + II - Ad + /(d)||. (3.23) 

Hence, 

/•oo 

/ \\dtiT)\\dT < +00, (3.24) 

J to 

11 



to 



which easily implies that as t — +00, d{x,t) converges in L^(fi). This and (13.13P indicate 
that 

lim \\d{t)-d^\\ =0. (3.25) 
Since d{t) is uniformly bounded in H^{Q) (cf. (13 ■31) ). by interpolation we have 

lim \\d{t)-d^\\Hi = 0. (3.26) 

t— »+oo 

On the other hand, uniform bound of d in H'^{Q) implies the weak convergence 

d{t) doo, in H^n). 

However, the decay property of the quantity A(t) (cf. Lemma [3^ could tell us more. 
Namely, we could get strong convergence of d in without using uniform estimates in 
higher order norm. To see this, we keep in mind that that 

\\Ad-Ad^\\ < \\Ad-Ad^-f{d) + f{d^)\\ + \\f{d)-f{do,)\\ 

< \\Ad- fid)\\ + \\f'iO\\L4d-d^\\L^ 

< \\Ad- f{d)\\ + C\\d-d^\\Hi. (3.27) 
The above estimate together with f l3.4p and (13.261) yields 

lim \\d{t) - dooWm = 0. (3.28) 

To finish the proof, we will show that Proposition 13.21 always holds true for the global 
solution d{t) to system (ll.ll) - (ll.5l) . Define 

tn = sup{t>tn\ \\d{-,s)-d^\\Hi < /5, V s G [t„,t]}. (3.29) 

It follows from (13.131) that for any e G (0,/3), there exists an integer N such that when 

n>N, 

\\di-,Q -d^Wni <€, (3.30) 
^{Eitn)-Eid^)f <e. (3.31) 

On the other hand, we can easily see that the orbit of d is continuous in H^. This is 
because we already know from ( 13.31) that d G L°°{0, +oc; H'^{Q)). As a consequence, 
d G L'^{t,t + 1;H'^{Q)) for any t > 0. The basic energy law and (13.23P imply dt G 
L2(t,t + l;L2(fi)). Thus, d G C{[t,t + l\;H\n)), for any t > (cf. [3]). The continuity 
of the orbit of d in and (I3.30p yield that 

tn > tn, for all n > N. 
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Then there are two possibilities: 

(i) . If there exists Uq > N such that ino = +00, then from the previous discussions in 
Case 1 and Case 2, the theorem is proved. 

(ii) . Otherwise, for all n> N, we have t„ < t„ < +00, and for all t G [tn,^n], E{doo) < 
£{t). Then from (I3.2ip with being replaced by and t being replaced by t„, we get 
from fIXaTD that 

/ {\\Vv{r)\\ + \\-M{r) + f\d{r))\\)dr<e. (3.32) 
Thus, it follows that (cf. (13:231 )) 

-c?oc|| < IM(tn) -rfooll + / IMt(r)|Mr 

< - rfooll + C / (||Vi;(r)|| + II - Arf(r) + /(rf(r))||)rfr 

< Ce, (3.33) 

which implies that lim„_,_|_oo ||'i(in) — c^ooll = 0. Since d{t) is relatively compact in H^iVt)^ 
there exists a subsequence of {d{tn)}^ still denoted by {d{tn)} converging to d^o in iJ^(fi), 
i.e., when n is sufficiently large, 

\\d{tn) - dooWm < P 

which contradicts the definition of t„ that \\d{-,tn) — doo\\m = P- 

Summing up, we have considered all the possible cases and the conclusion is 
proved. 



3.2 Convergence Rate 

In this part, we shall show the estimate on convergence rate (11.131) . This can be achieved 
in several steps. 

Step 1. As has been shown in the literature (cf. for instance, [9,30]), an estimate on 
the convergence rate in certain lower order norm could be obtained directly from the 
Lojasiewicz-Simon approach. From Lemma [2. II and (13.201) . we have 

j^{S{t) - E{d^)) + C,{S{t) - E{d^)y^'-'^ < 0, V t > to, (3.34) 
which implies 

£{t) - E{d^) <C{l + t)-^ yt>to. (3.35) 
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Integrating (I3.20p on (t, oc), where t > to, it follows from (13.231) that 

/OO POO 
\\dt\\dT<j^ {C\\\/v\\ + \\- Af + f{d)\\)dT <C{l+t)-^. (3.36) 

By adjusting the constant C properly, we obtain 

Mt) - d^W <C{l+t)-^, t>0. (3.37) 

Step 2. In Step 1, we only obtain the convergence rate of d (in L^). Unlike for the 
temperature variable in some phase-field systems (cf. [5,29] and references cited therein), 
although we have got some decay information for the velocity field v such that 

/OO 
\\\/v\\dT <C{l + t)-^, (3.38) 

it is not easy to prove convergence rate of v directly. This is because now v satisfies a 
Navier-Stokes type equation, which is much more complicated than the heat equation 
for the temperature variable in phase-field systems. As a result, one cannot easily obtain 
relation between || Vf || and Vt (in certain possible norm) from the equation itself. However, 
it is possible to achieve our goal by using the idea in [29], where we use higher order energy 
estimates and construct proper differential inequalities (cf. also [5,6,28]). Besides, in this 
way the convergence rate of d in higher order norm can be proved simultaneously. 

The steady state solution corresponding to problem p.ll) - p.5h satisfies the following 
system (cf. [15]) 

Voo-Vvoo-iyAvoo + VPoo = -V- (Vdoo© Vdoo), (3.39) 
V-v^ = 0, (3.40) 
^^oo ■ Vc/oo = Adoo- f{doo), (3.41) 
Vooir = 0, doo\v = dQ{x). (3.42) 

Lemma [T2] implies that the limiting point of system (|l.ip - (ll.5p has the form (0, d^o) € 
S. As a result, system ( 13.39p -( l3.42p can be reduced to 



|V4 



|2 



VPoo = -Vrfoo ■ Arfoo - V y^-^^j , (3.43) 
-Ad^ + f{d^) = 0, (3.44) 
doo\r = do{x), (3.45) 



where in (I3.43P we have used the fact that 



IV4 



|2 



V ■ (Vrfoo Vrfoo) = V (^^^ ) + Vrfoo ■ Ad^. 
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Subtracting the stationary problem (I3.43p - (l3.45p from the evolution problem (ll.ip - dl.Sp . 

we get 

Yivrf|2\ /|V4 "~ 



"oo I 



Vt + v-Vv -uAv + V{P - Po., ■ • , , 2 

= -Vrf ■ Arf + Vdoo ■ Adoo, (3.46) 

\/-v = 0, (3.47) 

dt + vVd = A{d- rfoo) - fid) + f{d^), (3.48) 

(d-doo)|r = 0. (3.49) 

Multiplying I KM by f and (IXiSll by -Ad + f{d) = -A{d - d^c) + fid) - f{doo) respec- 
tively, integrating on f2, and adding the results together, we obtain 

4 filkf + il|Vrf-Vrfooir+ / F{d)-F{doo)-f{do.){d-d^)dx 



dt V2' 2' 

+u\\Vvf + \\Ad-f{d)f 



n 



= {v, Vrfoo ■ Adoo) 

= (t;, Vc/oo ■ (Ac/oo - f{doo))) + {v ■ Vdoo, -/(c^oo)) 

= 0. (3.50) 
Multiplying fl3.48l) hy d — d^o and integrating in fl, we have 

~\\d~d^f+\\V{d~d^)f = -{vVd,d-d^)-{f{d)-f{d^),d-d^) :=/i. (3.51) 
The right hand side can be estimated as follows 

|/i| < \\v\\L4^d\\L4d-d^\\ + \\f'{mL4d-doo\\l^ 



ooll)^ 



< C\\Vv\\\\d-d^\\+C{\\V{d-d^)p\\d-d^p + \\d-d, 

< ei\\Vvf + ^\\V{d-doo)f + C\\d-doof. (3.52) 

Multiplying (13.511) by a > and adding the resultant to (13.501) . using (I3.52p we get 

d (\ 



c/t V 2 '' " 2 



l\\Vd-Vdo^f + ^\\d-d^f+ [ F{d)dx- [ F{d^)dx 
2 ^ Jn Jq 



n 

2 



f{d^){d - d^)dx ]+{v- aei) \\Vvf + \\Ad - f{d)f + ||| V(rf - d 



2 , "nw.. . ^||2 



< Ca\\d-doo\\. (3.53) 
On the other hand, by the Taylor's expansion, we have 

F{d) = F{d^) + f{d^){d - d^) + /'(0(rf - d^)\ (3.54) 
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where ^ = ad + {1 — a)doo with a G [0, 1]. 
Then we deduce that 



F{d)dx — / F{doo)dx + / f{doo)doodx — / f{doo)d dx 
n Jn Jn Jn 



f'iOid-d^fdx 



< ii/'(oiUociid-doof <c2iid-dooir. 

Let us define now, for t > 0, 

yit) = ^l|t;(t)f + i||Vrf(t)-Vrfoof + |lM(t)-^^ooir+ /" 

- / f{doo){d{t) - doo)dx. 
Jn 

In ( 13.53P and (I3.56p . we choose 



(3.55) 



F{doo)dx 
(3.56) 



As a result, 



a > 1 + 2C2 > 0, ei = —. 

4a 



y{t) + - rfooll' > -dl^^ir + \\d - rfooll^i). 



Furthermore, we infer from (13.571) that for certain constants C3, C4 > 0, 



j^y{t) + Csvit) < C4d - doof < c{i + tyT 



29 
-29 



As in [28,29], we have 



y{t) <C(l + t)-^, Vt>0, 



which together with ( 13.57p implies that 



\v{t)\\ + \\d{t) - d^Wni < C(l + 1)-—, V t > 0. 



(3.57) 



(3.58) 



(3.59) 



(3.60) 



Step 3. In the last step, we proceed prove the convergence rate in higher order norm. 
In Section 3.1, it has been proven that, once we could obtain the uniform bound of d in 
if^, we are able to obtain strong convergence of d in instead of weak convergence. 
By reinvestigating the higher order energy estimate for the subtracted system ( 13.46^ - 
f l3.48p (cf. also Lemma [3TT]) . we can obtain a further result, which provides the same rate 
estimate of {v, d) in x H"^ as (ITBnD . 

In what follows, we just perform the estimates for classical solutions. Taking the time 
derivative of A{t), we obtain by a direct calculation 

1 d 



2dt 



A{t) + {\\Svr+\\V{Ad-f{dW 
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(Sv^vVv) - if'id)iAd~ f{d)),Ad- f{d)) + 2 [ {Ad-f{d)\,vi^d^Jx 

- j Vn -Vd^Ad- f{d))dx 

{Sv, vVv)- {nd){Ad - f{d)), Ad - fid)) 

+2 / (Ad- f{d)),^vi^id-d^),Jx-2 [ {Ad- f{d))WvW^d^dx 
Jn Jn 

- I -VdiAd- f{d))dx 



h + h + h + h + Jg. 



(3.61) 



In the above, we use the fact that {Sv,Vt) = {—Av,Vt), which follows from Vt G H. 
Noticing that we have got uniform bounds for ||f and ||(i||//2 before (see Lemma [X2|) . 
in what follows we estimate (i = 2, 6) term by term. 



I/2I < \\Sv\\\\v\\l4^v\\l^ < C\\Sv\\{\\Vvp\\vp)i\\Avp\\Vv\ 
< C\\Sv\\\\Av\\-^\\vf^ < e2\\Svf + C\\vf. 



(3.62) 



Since 

< 
< 
< 
< 

we have 



||V(Arf- Adoo)|| 

\\V{Ad-f{d))\\ + \\V{f{d)-f{d^))\\ 

||V(Arf - fid))\\ + \\f'id)iVd - Vd^)\\ + Uf'id) - /'(rfoo))Vrfoo|| 

||V(Arf-/(d))|| + ||/'(rf)|U^||(Vd-Vdoo)|| + ||/"(OI|Loo||d-cioo|U4||Vrfoo||L4 

||V(Arf-/(d))||+C||rf-rfoobi, (3.63) 



I/3I < \\r{d)\\L^\\Ad-f{d)f<C{\\Ad-Ad^\\ 

< C\\Ad - Arfoof + CWrml^Wd - d^f 

< C\\V{Ad-Ad^)\\l\\d-doo\\-^+C\\d-do 

< e2\\V{Ad~fid))f + C\\d-doof. 



\\fid)-f{do 



(3.64) 



Next, 



I/4I < ||V(Arf-/(rf))||||Vt;|U4||V(d-doo)||L4 

< e2\\ViAd-fid))r 

+C(||A^;||||Vi;|| + ||Vi;f)(||A(rf-rfo.)||||V(d-rfo,)|| + ||V(rf-rfoo)f) 

< £211 V(Arf - f{d))f + e^WSvf + C\\A{d -d^)f + c(l + ^] \\V{d - d^)^ 
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< 2e2\\ViAd- fid))f + e2\\Svf + C\\d-d^\\l^. (3.65) 



\h\ < \\Ad- f{d)\\L4^v\\L4^^d^\\ 

< C\\V{Ad ~ f{d))\\{\\Av\\^\v\\-^ + \\v\\) 

< £2||V(Arf - /(d))f + -i\\Av\\l\\vf^ + \\vf) 

< e,\\V{Ad-f{d))f + 62\\Svf + c(\ + -] \\vf. (3.66) 



\h\ < ||V7r||||Vrf|U4||Arf-/(rf)|U4 

< C\\Sv\\\\ViAd- fid))\\-^\\Ad- fid)f^\\d\\H2 

< e2\\V{Ad - /(rf))f + e2\\Svr + ^\\Ad - f{d)f 

< 2e2\\V{Ad-f{d))f + e2\\Svf + ^\\d-d^\\l^. (3.67) 
Taking 62 sufficiently small, we deduce from fl3.6ip - fl3.67p that 

j^A{t) + {\\Sv\\' + II V(Arf - /(rf))f ) < C(||t;f + ||d - d^Wl,). (3.68) 

Using the Poincare inequality for Ad — f{d) whose trace on F is and Lemma 12^21 we 
can conclude from f l3.68p and fl3.60p that 

—A{t) + CA{t)<C{\\vf + \\d-doo\\m)<C{l+t)-^, Vt>0. (3.69) 

Again, by the Gronwall inequality, we have 

A{t) <C{l + ty^, Vt>0, (3.70) 

which yields 

\\Vv{t)\\ + \\Ad{t)- f{d{t))\\<cii + ty^, yt>o. (3.71) 

Recalling fJS^ZD, it follows from fl3?7Tll that 

\\Ad{t) - Ad^W <C{l + t)-^, Vt>0. (3.72) 

Summing up, from f l3.60pfl3.7ipfl3.72p we can deduce the required estimate p.l3p . The 
proof of Theorem 11.11 is complete. 
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4 Results for Three Dimensional Case 



The results proved in previous section hold true for global classical solutions to system 
f ll.ip - (ll.5l) in 3-D case. In what follows, we show the convergence to equilibrium for two 
subcases considered in [15] (ref. Theorem B and Theorem C therein) that existence of 
global classical solution was proven. In particular, we answer the question of uniqueness 
of asymptotic limit ofd (cf. [15, Remark, page 32]) and provide a uniform convergence rate. 

Case I: Initial Data Near Absolute Minimizer of E. 

The following result has been proven in [15, Proposition 5.2]. 

Proposition 4.1. There is an Eq G (0,1) depending only on z/, A,7,f2 and f with the 
following property: Whenever 

u\\VvfiO) + X^\\Ad-f{d)f{0)<eo, 

either 

(1) Problem p.ip - p.5l) has a unique classical solution {v,d) in Q x (0, +oo) 
or 

(2) there is aT^ E (0, +oo) such that 

E{T,) < E{0) - Eo, 

where 

E{t) = \\vf + X\\Vdf + 2\ I F{d)dx. 

Jn 

Moreover, in case (1), one has 

\\v{t)\\H^n)-^Q, ||Arf-/(rf)|| ^0, ast-^+oo. (4.1) 

Before proving the convergence result corresponding to Theorem 11.11 we turn to the sec- 
ond case. Later we shall prove our result in a unified way. 

Case II: Arbitrary Initial Data with Large Viscosity. 

It has been proven that for any initial data Vq G H^{Vt)^ do G if^(fi), if the viscosity u 
is "large enough" (see below), problem ( ll.ip -( fT75l l admits a unique global classical solution 
(cf. [15, Theorem B]). As pointed out in [15], when the dimension is three, the size of 
viscosity z/ plays a rather crucial role while the other constants A, 7 do not, as long as A, 7 
are positive constants. Thus we shall assume A = 7 = 1 for the sake of simplicity. The 
following high order energy estimate can be obtained (cf. [15, (4.13)]). 
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Lemma 4.1. In the 3-D case, the following inequality holds for classical solution {v,d) 
to problem (frTD-ffLBD 

~A{t) < - {u - Ku'^A) \\Avr-(^l - ^^iVi^^^ Vt>0, (4.2) 

where A = A + 1 = ||Vf |p + \\Ad — /(c?)p + 1 (cf. (I3.ip ) and K is a positive constant 
depending on /, \\vo\\, WdoWn^n)- 

When the viscosity u is assumed to be properly large, based on the above lemma, we 
can not only show that the global solution {v,d) is uniformly bounded (as in [15]) but 
also the quantity A{t) decays to zero in time. 

It follows from (13.61) that 



A{T)dT < A{T)dT + 1 < M, V t > 0, (4.3) 

where M > is a constant depending only on ||fo||, IMoIIhi- Then we have 
Lemma 4.2. // 

>K (^A{0) + 2KM + 4:M^ (4.4) 

then the unique global solution to problem 01.11) - (1 1.51) satisfies the following uniform esti- 
mate 

\\v{t)\\Hi + \\dmH^<C, Vt>0, (4.5) 
where C is a constant depending on f,^, \\vo\\H'^{n)i \\do\\H^(n)- Furthermore, 

lim {\\v{t)\\H^ + II - Ad{t) + f{dm) = 0. (4.6) 

t— >+oo 

Proof. Proof of existence and uniqueness of the global solution has been given in [15]. 
Next, we show the uniform bound (14.51) . Take u large enough that (14. 4p is satisfied. Then 
by (|4.2p . there must be some Tq > such that 



u-Ku-^A{t)>0, i-:^M>o, 



for all t G [0,To]. Moreover, on [0,To], 

|i(t) < 2KA{t). (4.7) 

Denote = supTg. First we show that T* > 1 by a contradiction argument. 
If < 1, then 

Mt*) < ^(0) + 2K I A{t)dt < i(0) + 2KM. 
Jo 
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On the other hand, from the definition of T*, we have 



V < max{KA{T,), K'A'{n)} < K{A{0) + 2KM) + K\A{<S) + 2KM)\ 

which contradict (14 ■4p . 

Next, if T* < +00, f l4.3p implies that there is a ti G [T* — |,T*] such that 

A{ti) < AM. (4.8) 

As a result, 

A{T,)<AM + 2K I A{t)dt<AM + 2KM. (4.9) 
Jtl 

Again from the definition of T*, we have 

z/ < ma.x{KA{T,),K^A^{T,)}, 

which together with (14.91) yields a contradiction with (14.41) . 
Therefore, for all t > 0, (lO) holds. Namely, 

j^A{t) < 2KA{t) +2K < KA\t) + 3K. (4.10) 

Due to (13. 6p . we can conclude l \AM and ( 14.61) following the similar argument in the proof 
of Lemma 13. 2[ □ 

Remark 4.1. Generally speaking, (14.41) only provides a sujficient condition on the large- 
ness of viscosity v, which ensures the existence of global solution to problem (ll.ip - (ll.5l) . 
It may not be an optimal lower bound for all possible u. 

Based on above results, now for both cases I and II, one can argue exactly as in Section 
3.1 to conclude 

lim {\\v{t)\\m + \\d{t) - dooIkO = 0. (4.11) 

Then we are able to proceed to show the estimate on convergence rate for both two 
cases. To this aim, we check the argument for 2-D case step by step. By applying 
corresponding Sobolev embedding Theorems in 3-D, we can see that all calculations in 
Section 3.2 are valid for our current case (with minor modifications). Hence the details 
are omitted. 

We complete the proof for Theorem 11.21 and Theorem 11.31 
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5 Further Remarks 



We remark that our approach used in this paper are valid for some other model systems 
for nematic liquid crystal flows in the literature and similar convergence result can be 
proved. 

(1) A model with changing density 

Recently the following problem was considered in [20]. 

Pt + S/-{pv) = 0, p>0, (5.1) 

{(w)t + V ■ {pvQv)-u^v + VP = -XV ■ (VdQVd), (5.2) 

\/-v = 0, (5.3) 

d, + vVd = 7(Ac/-/(rf)), (5.4) 

in X (0, oo), where Q C ]R"(n = 2,3) is a bounded domain with smooth boundary T. 
p{x, t) is a scalar function denoting the density of the fluid. The above density-dependent 
liquid crystal model is subject to the following initial condition 

p\t=o = Po{x) > 0, ipv)\t=o = qo{x), d\t=o = do{x), ior x e fl, (5.5) 

and the boundary conditions: 

v{x,t) = 0, d{x,t) = do{x), for (x,t) G r X M+. (5.6) 

Problem (15.1l) - (l5.6p can be viewed as a generalization of our problem f ll.ip -( fT75l ). It 
enjoys some important properties as for p.ll) - p.5p . In particular, we have the following 

basic energy law (see [20]) 

Jt i [1^^''^' + + ^ / ^^^)^'^) =~i + ^^'^^ ~ ^^^^''^ ^^-^^ 

In [20], the authors proved the existence of the weak solution to incompressible liq- 
uid crystal system ( 15.1l) - (15.6p under certain compatibility condition on the initial data. 
There they considered the general case for the density, namely they only required the 
initial density to be nonnegative. As a result, one can only expect the density to be 
nonnegative for all time and vacuum state may occur. In this case, it is very difficult to 
prove corresponding results to [15, Theorem B, Theorem C] where the global existence 
and uniqueness as well as asymptotic behavior of classical solutions were obtained. 

However, if we assume in addition that the initial density is a bounded positive func- 
tion, i.e., there are two positive constants p and p such that 

< p < po{x) <p, ^ X eQ. (5.8) 
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Then by virtue of the comparison principle (cf. [20]), we have 



< p < p{x,t) < p, Vt>0. 



(5.9) 



In this special case, we can check that under suitable assumptions on the initial data 
(for instance, vq G Hq, do G i/^), parallel results to [15, Theorem B, Theorem C] can be 
achieved. Besides the basic energy law (15 .Tp . due to uniform upper and lower bound (15. 9p . 
one can proceed to get proper high order energy law similar to Lemma [3T] (2-D case) as 
well as Lemma HT] (3-D case). For instance, denote 



we can show that 

Lemma 5.1. In the 2-D case, the following inequality holds for the classical solution {v, d) 



-Mt) + m\Svf + II V(Ad - /(d))||2) < KM^t) + 1), Vt > 0, (5.11) 



where Ki,K2 are constants depending on f,^l, ||fo||, \\do\\m{Q),i^, ^,1, P, P- 

Corresponding results in 3-D case (cf. Lemma [HI]) can also be obtained. The proofs 
for these results follow from the same sprit of those in [15] with some proper modifications. 
Hence, the details are omitted here. 

Remark 5.1. Different from system (ll.ip -( fT75ll . where the density is assumed to be a 
constant, in order to get the high order energy law, we deal with the time derivative of 
a modified quantity A{t) with weight p instead of A(t). This is due to the mathematical 
structure of (I5.ip - (15.6p . in which the density variable is involved. Because of the uniform 
upper and lower bounds of the density (15.91) . one can check that A{t) plays a similar role 
as A{t) for system (fTTD-fra. 

Based on the facts obtained above, we are able to prove the corresponding convergence 
results (cf. Theorem 1.1-Theorem 1.3) for system (I5.ip - (l5.6p . following the argument in 
the previous sections. We leave the details to the interested readers. 

(2) A model with free-slip boundary condition 



A{t)=p\\vvr+x\\Ad-fid)r, 



(5.10) 



to problem (15. II) - (15.61) 



vt + v - Vv - udivD{v) + VP 



AV ■ ( Vrf Vrf) 



(5.12) 
(5.13) 
(5.14) 



V ■ V 







dt + vVd 



7(Arf-/(rf)) 
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in i7 X (0,00), where Q C W^{n = 2,3) is a bounded polygonal domain (with piecewise 
smooth boundary). D{v) = j(Vv + (Vv)'^) is the stretching tensor. We consider the 
system (I5.12p -( l5.14p subject to the initial conditions 

v\t=o = Vo{x) with V ■ fo = 0, d\t=o = do{x), for x G fi, (5.15) 

and the free-slip boundary conditions: 

v-n = 0, (V X f ) X n = 0, d^d = 0, for {x, t) eV x R+, (5.16) 

where n is the unit outer normal vector to the boundary T. 

As has been pointed out in the recent paper [18], the free-slip boundary condition 
(I5.16P indicates that in the liquid crystal flows, there is no contribution from the director 
field d to the surface forces. Boundary condition (15.16p seems to be more appropriate for 
some types of flow in the bulk of a liquid crystal configuration. On the other hand, it 
allows people to construct more efficient numerical schemes for the numerical simulations 
for liquid crystal flows (cf. [18]). Comparing with system p.ip -f fT75ll . the influences of the 
corner singularities is less severe with free-slip and Neumann boundary conditions than 
the Dirichlet boundary conditions. 

Basic theoretical analysis on problem (I5.12p - (l5.16p has been done in [18], where the 
authors proved global existences of weak solutions as well as regularities and global ex- 
istence/uniqueness of classical solutions. In particular, although the boundary condition 
(I5.16P plays a significantly different role in the calculation, proper high order energy law 
similar to Lemma [SH] could still be obtained (cf. [18, Lemma 4.1]). The same convergence 
results for system p.ip - (11.5p obtained in the present paper can be shown true for problem 
fl5.12p -( !5.16p . by adapting the argument here. We thus omit the details. 
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